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We chart the theoretical basis of radiative decays of the Higgs boson, H — > 77 and Z'y, for measuring 
the spin of the Higgs particle. These decay channels are complementary to other rare modes such as 
real/virtual Z-boson pair-decays. In systematic helicity analyses the angular distribution for zero- 
spin is confronted with hypothetical spin-2 ± and higher assignments to quantify the sensitivity. 

1. After the discovery of the Higgs particle, the properties must be examined experimentally to identify the 
particle as the element proper of the Higgs mechanism for breaking the electrowcak symmetries [l| [for recent 
general reviews see Refs. Q]. The dynamical steps include the confirmation of the scalar spin-zero character of 
the particle 1 . Higgs-boson decays to pairs of Z-bosons 2 , cf. Refs. [3-5] and additional references listed there, 
are promising candidates for examining the Higgs spin at the LHC, supplemented by other measurements in the 
Higgs-strahlung processes at LHC Q and e + e~ linear colliders Q. Exploiting the primary Higgs-boson search 
channels, 77 and WW decays provide another self-evident tool for spin measurements 0, 0|. [For a variety of 
methods see the recent literature in Refs. 10]]. All these channels are particularly difficult to analyze for low masses 



of the Higgs boson where the 6-pair decays are overwhelming. However, with about 126 GeV this is precisely the 
mass range, where a new boson has been discovered by the LHC experiments 



11 



13 1, consistent with analyses 



of electrowcak precision data |14| and compatible with Higgs patterns within the framework of the Standard Model 



In this study we examine to which extent radiative decays 16l-ll8j 



X 



H — > 77 and Z7 



(1.1) 



for masses of the order of 126 GeV 
generates a sample of order 10 4 events in 77 decays for an LHC luminosity of 



can be exploited to determine the spin of the Higgs boson. Even though these decays are rare, with BR ~ 2 • 10 3 

the large number of light Higgs bosons produced in gluon fusion 

100 fb _1 , and a similar number 
of ^7-decay events [though reduced finally by the branching ratio for leptonic e, fj, decays of the Z-boson]. These 
channels nicely supplement the other rare ZZ* channels. In the present report we perform general helicity analyses 
of the radiative decays to quantify the sensitivity to zero-spin of the Higgs boson in angular distributions. 



2. The zero-spin character of the Higgs boson reflects itself in the isotropic decay distribution of the radiative decays 
in the rest frame, i.e. 

1 <ff 77 , Z7 _ i ^ i/2 



g?cos 9 



with denoting the polar angle of the 77 and Zj axis, defined, for example, in the Higgs-boson rest frame with 
regard to the LHC beam axis, Fig.[lja), and with the normalization including the proper statistical factor of the 
two-photon state. The axis can be reconstructed experimentally in the exclusive 77 and Z7 final states. Since the 



1 This letter expands part of a talk delivered by S.Y. Choi at the APCTP 2012 LHC Physics Workshop at Korea, Seoul, Korea, August 
7-9, 2012. 

2 In continuum final states of Higgs-boson decays H — > Z + [//] the fermion pair [//] is not necessarily decay product of a virtual Z* 
state so that general experimental analyses are not constrained by Bosc symmetry. 
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Z-boson is produced in helicity ±1 states, the U distribution in leptonic Z decays is of the familiar 1 + cos 2 9i form, 
6i being the polar lepton angle in the rest frame of the Z-boson. 
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FIG. 1: (a) Kinematics of radiative Higgs-boson decays in gluon fusion; (b-d) Angular distributions 0/77 and Z7 axes in the 
rest frame of the subprocesses: the flat Higgs signal compared with the distributions of the background events with the invariant 
qq energy set to Mh = 126 GeV in (b) and all potential spin-2 distributions corresponding to the spin component along the 
collider axis in the rest frame of the Higgs boson, and the difference of the helicities of the vector bosons in the decay final states 
[all distributions normalized over the interval | cos 9| < 1/V2/. The upper indices refer to the allowed parity associated with the 
distributions in 77 decays. Allowed parities of Z7 states are denoted in round brackets either if unique for Z'y or if different 
from 77 final states. [Note that Z7 states do not conform with the gg initial states in contrast to 77 states.] 



The flat distribution (|1.2[) is unique to zero-spin of the Higgs boson. For any other spin assignment J the 
distribution would be given by the Wigner functions ~ l^m.A -A' (®)| 2 an d l^m,A z -A' (®)l 2 > wnere m denotes the 
S z spin component being either or ±2 for Higgs-boson production in gluon fusion, while A 7 — = 0, ±2 and 
Xz — A' = 0, ±1, ±2. The explicit form of the Wigner d-functions d^ A (0) may be read off the tables in Ref. [23). 



The general parity-invariant polar angular distribution can be cast into a compact form for any Higgs-spin J decays 
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to 77 and -Z7 in gluon fusion: 



(1.3) 



and 



1 d 2 <r[Zj] 
a dcos <d dcost 



3(2J + 1) 
16 



x J y' J v£ + x J y' J 2 v J 02 + x 2 J y'iv( + x 2 J y J 2 v{ 2 (i + cos 



'JM> J r<J 



(1.4) 



respectively, with the squared Wigner functions 2>^ A = j J2[dm\(®)] 2 symmetrized by the summing over all signs 



±m, ±A of the spin/helicity components [index counting twice]. The functions read explicitly for J = and J = 2: 



L \)o — 1 



V 2 



(3 cos 2 6 - l) 2 /4 



V 2 20 =V 2 )2 = 3 sin 4 0/8 



V 2 



V 2 10 = 3 sin 2 6 cos 2 6/2 



V 2 22 = (cos 4 6 + 6 cos 2 + 1)/16 



V 2 



10 

"D 2 



sin 2 6(1 + cos 2 6)/4, 



(1.5) 



the functional forms displayed in Fig.QJc,d). The [non- negative] reduced production and decay helicity probabilities 
X for gluon fusion, y for 77 and y' for Z"f final states, are mo del- dependent parameters, obeying the sum rules 



Xn 



■ xi = 1 



yi + yi = 1 and y'J + y[ J + y' 2 J = 1. 



(1.6) 
(1.7) 



[The formalism can easily be transferred to qq production with S z = or ±1 by substituting Xf for X 2 and T>^ k for 
T> 2k (fc = 0, 1, 2) in the cross sections and in the X sum rule. Note however that the rate for signal Higgs production 
in qq collisions is negligibly small for light quark beams. The initial states gg and qq mix incoherently in the most 
general configurations.] 



The polar angle of the Z decay distribution can easily be integrated out, equivalent to the substitutions 
(1 + cos 2 6» £ ) -> 8/3 and sin 2 t -> 4/3. Since the reduced helicity probabilities are non-negative, the coefficients of 
both the (1 + cos 2 61) term and the sin 2 Bg term inevitably generate non- vanishing maximum/ minimum ±cos 2J 6 
terms. Thus, observing [in the experimentally ideal case] that the angular distribution is independent of cos 6, proves 
unambiguously the spin-zero character of the Higgs particle. At the same time the sin 2 9g term in Zj final states is 
predicted to be absent, providing an independent cross-check. 
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TABLE I: Selection rules for Higgs parity following from observing the polar angular distribution of a spin-J Higgs state in the 
process gg — > H — > 77. 



The observation of spin states in gg — >• H —> 77 allows [partial] conclusions also on the parity of the states. 
From Bose symmetry and parity symmetry of the helicity amplitudes T\ = (— ^-) J T\ \! = ^C - 1) J T-\ 1 -a > 
referring separately to initial and final states the selection rules presented in Tab. Q] can easily be derived [see 
also Ref. [241]] . complementing global rules noted earlier in the literature. Scalar and pseudoscalar Higgs bosons 
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cannot be discriminated in the 77 decay mode, neither even/odd parity by observing T>q Q , and spin correlation 
effects 0, 0, [3] must be exploited for discrimination. But observing any state with helicity difference AA = 2 
in initial or final state determines unambiguously the even-parity character of the Higgs boson. Analogous rules 
apply also to Z7 final states in gluon fusion and 77 final states in qq annihilation. Either the first or the second 
index AA in the V functions coming with the production or decay amplitude, respectively, is restricted in paral- 
lel to the rules in Tab.Uwhile the respective companion index is unrestricted apart from the standard spin constraints. 

The 77 channel is described by two- by- two independent probabilities for production and decay, Xq 2 and y^ 2 - 
Popular choices for experimental simulations are the complementary {J; 00} 'scalar-type' and the {J; 22} 'tensor-type 
assignments ': 

'scalar-type assignment' : Xrf = yg = 1 and X 2 J = y£ = [J > 0] (1.8) 
'tensor-type assignment' : X^ = y^ = and X( = y 2 J = 1 [J > 2] , (1.9) 

supplemented by the 

'mixed-type assignment' : Xq = y 2 = and X 2 = 3^ = 1 [J > 2] and '1,0' interchanged, (1.10) 

for the signal J = and the hypothetical alternatives J > 2. [The 77 coupling of the tensor- type assignment is 
equivalent to the KK graviton coupling in d = 5 scenarios [sj].] The configurations can be exploited in two ways: (i) 
One of the two scalar- or tensor-type configurations for J > 2 for instance, is sufficient to prove that the spin-zero test 
of the signal is non-trivial; (ii) However, to prove experimentally that the spin-2 assignment is not realized, the three 
configurations, which are mutually independent, must necessarily be shown absent. Not observing the double index 
{J; 00}, the state { J~} is ruled out, while neither observing distributions carrying at least one index 2, the state { J + } 
is ruled out, too. Thus any spin J > 2 can be rejected for both parities ± by angular analyses. - Due to potential lon- 
gitudinal Z polarization accounted for by y'f, the Z7 final state is described by three independent decay probabilities. 



Disregarding J = 1 , as forbidden by the Landau- Yang theorem in 77 decays [24j, |2jj , we will choose J = 2 for 
illustration. [Part of the distributions have also been noted in Rcfs. None of the possible helicity states would 

generate a flat distribution like spin-zero Higgs-boson decays. It is shown in Fig.Q] how the flat signal distribution 
contrasts with the distributions generated by hypothetical spin J = 2 assignments of both even and odd parity, 
and the angular distribution of the backgrounds as well [to be analyzed next]. In the final section we will compare 
the spin-0 distribution, assigned J = and m = A 7 ,z — A^ = 0, specifically with the scalar-type {2; 00} and the 
tensor-type {2; 22} distributions, representing the state J = 2 with the two parities ±. The {2; 00} scalar-type 
distribution is pronounced in the center like spin-0 after angular cut, and rises in the forward/backward directions 
like the continuum backgrounds, thus providing a non-trivial analogue to be discriminated experimentally from the 
Higgs signal. The {2; 22} tensor-type distribution rises monotonically to the left and to the right of the center, 
providing also a valuable discriminant. Both types must be ruled out necessarily to reject experimentally the spin-2 
assignment for even and odd parity. 



A first global comparison between spin-0 and all the spin-2 distributions is offered by the moments of the 
polar-angle distributions of the 77 and Z7 event axes 3 noted in Tab.|TTJ The first moments of the spin-2 distribu- 
tions are characteristically different 4 from the spin-0 distribution and may provide early information on the Higgs spin. 

3. The large continuum background generated in pp(qq) — > 77 and Zj processes [and, to a lesser extent, loop-induced 
gluon fusion] requires stringent cuts in order not to dwarf the signals. The angular characteristics allow to reduce 



3 Angular asymmetries [i| can equivalently be used in global analyses. 

4 The assignments {2; 01} and {2; 21} can also be discriminated from {0; 00} by identifying the Z-dccay angular distributions sin 2 t 
versus (1 + cos 2 9g). 
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TABLE II: Ratios of the first over zeroth moments of the polar angular distributions of the 77 and Z-y axes, | cos0|, for zero 
Higgs spin and for spin 2 of a hypothetical resonance, 'cut ' denotes the theoretical cut \cosQ\ < l/\/2 on the polar angle of the 
event axes. 



the continuum backgrounds considerably. 



The angular distribution of the background events is strongly peaked in the forward and backward directions as a 
result of the t and it-channel exchange mechanisms, in contrast to the flat distribution of the signal. In the rest frame 
of the parton system, cutting out the singular forward and backward directions — s- and n [so long as the cross 
sections are not regularized properly]: 



(1(7 



dcosO 

da 
dcos 9 



[qq 77] 



[qq Zj] = 



2na 2 
3s 

2na 2 
3s 



1 



sin 2 6 



[1 



■e] 



+ a 2 q ] [l-MZ/s] 



1 



sin 2 e 



1 + cos 2 e 



4M|s 



Ml? 



(1.11) 



(1.12) 



for the invariant parton energy y/s = Mu, electric and weak quark charges denoted by Q,v,a |29|, |3fJ; the helicity 
decomposition is familiar from electron-positron collisions, cf. appendix in Ref. [3lj |. The parton subprocesses are 
peaked at small angles, regularized by the strong interaction scale Aqcd- For the given Higgs mass, the angular 
distribution of gg — > 77 |32| is surprisingly close to the qq process, cf. Fig. [TJb). Induced by radiative return 33[ in 
j?7, and self-evident in 77, the 7's and the Z-bosons are traveling primarily along the LHC beam axis. By restricting 
I cosO| to less than cos0 ctl t < l/y/2, the signal is reduced only modestly, but the background strongly. 



4. For illustration of the J p sensitivity, we present a set of r oug h theoretical estimates, rescaled from available 



experimental data in 77 11[ or based on simulations in Zj [34| . One photon with transverse energy e porp in 



excess of 25 GeV was required for Z7 decays, and photons in excess of 40 GeV in 77 decays. Numerically, using 
cos 6 = [1 — iep^pMfj / (M^j — M|) 2 ]3 and [1 — 4ep 0rp /M^]5 for the corresponding polar angles in the Higgs rest 
frame, this is approximately equivalent to the restrictions | cos0| < 0.77 and 0.55, close to the theoretical cut \ 
used in the previous section. The Z-boson was assumed to decay leptonically to e, /i pairs. 



The error estimates for the 77 final states were performed by rescaling the background event numbers of Ref. 11 1 
in energy, using the theoretical energy dependence of qq —> 77 [ignoring mis-identification from jets in this simplified 
theoretical estimate], and by raising the luminosity to 100 fb" 1 , leaving us with 4.6fc [146fc] events after cuts for 
signal [background], in rough agreement, within a factor two with Ref. 35j, after inserting the proper _ftT-factor. The 
theoretical angular distribution of the photons in the center range left by the cuts was used as well according to 
Fig.[T] We have adopted the experimental efficiency of 40% and a resolution of ±2 GeV. Experimental refinements 
such as smearing effects etc. have not been considered in this coarse theoretical picture. 



In the same way as Ref. [34[ we determined the background event number from the cross section of the qq — > Z7 
process which dominates the background compared with gg collisions and the Drell-Yan process including final-state 
photon radiation. Similarly to the parameters in Ref. 3J] we have adopted the values of 3 GeV for the mass resolution 
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FIG. 2: Left: The angular distribution of the decay axis for spm-0 Higgs-boson decays H — > 77; the spin-0 expectation is 
compared with the scalar and tensor-type decay distributions of a hypothetical spin-2 particle [theoretical errors only for the 
LHC luminosity = 100 fb~ 1 , event numbers scaled up from fl~k l; Right: The same distributions for H — > Z-y decays [LH-LHC 
luminosity = 3 ab^ 1 ]. 



and 0.13 for the efficiency 5 . The production cross section of the signal was recalculated theoretically by including the 
large -ftf-value in Higgs-boson production pp — > H. Finally, 1.2k [51k] signal [background] events were predicted for a 
luminosity of 3 ab _1 . 

These theoretical estimates of rates and errors, not including detailed experimental refinements, should serve only 
as a rough illustration of theoretical expectations for Higgs spin analyses in radiative decays. 

In the framework defined above, the results for the signals and the size of the roughly expected backgrounds are 
shown in Fig. [2] for H — > 77 and H — > Zj in the theoretically cut cos0 range. The solid lines are the [ideal] signals 
with very high statistics, the error bars, based on the event numbers defined above, represent the theoretical estimates 
of background fluctuations \J~B, contaminating the angular distributions of the signals S [<B], These distributions 
are compared with the expected decay characteristics of a hypothetical spin-2 particle, even/odd parity, derived for 
the representative angular tensor-type and scalar-type distributions in Eqs. (|1.5[) . Note that the two distributions 
are mutually complementary to each other. The first |cos0| moments [normalized to zeroth moments] are in the 
ratio (0.52 ±0.06)/(0.27± 0.05) for the tensor/scalar assignment in 77 final states, and (0.36 ± 0.10)/(0. 18 ± 0.08) in 
Zj final states. Even if added up, the resulting flatfish behavior is fractured by the non-zero 2?20 = ^02 contributions. 

The small leptonic Z branching ratio together with the increased background cross section render experimental 
Zj analyses much more demanding than 77 analyses, and a large increase of luminosity is required. 

Studying experimentally the 77 and Zj processes outside the Higgs mass window will yield a good understanding 
of the background shapes and normalizations. At the expense of a y/2 increase of errors one could define a control 
region with a lower Higgs-typc mass window to determine the shapes and use Monte Carlos to extrapolate from the 
control region to the signal region. 



5 These parameters were extracted by means of PYTHIA [36| and ACERDET [37( . with support by D.Zcrwas gratefully acknowledged. 
The transverse momentum of the photon was chosen in excess of 25 GeV. In Z decays to leptons, electrons or muons were reconstructed 
with transverse energy /momentum of more than 25 GeV, separated from the photon by AR > 0.7, where AR is the square-root of the 
sum of the squares of the azimuthal angle and pseudo-rapidity differences. The invariant mass of the lepton pair was chosen compatible 
with the mass of the Z boson within 5 GeV. The reconstructed mass of the Higgs boson was allowed in a window of 3 GeV centered on its 
nominal mass. The second window can be chosen more restrictive as the Higgs boson width is negligible compared to the experimental 
resolution, whereas for the Z boson this is not the case. The resulting efficiency of 0.13 is close to the value reported in Ref. [34H when 
interpolating slightly different parameters. 
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5. Summary: Dynamical characteristics of the Higgs boson in the Standard Model are particularly difficult to 
analyze experimentally in the mass region around f 26 GeV since the overwhelming decays are b decays. In this report 
we have analyzed the theoretical potential of two decay modes, H — > 77 decays and H — > Z7 decays [the latter 
statistically more remote], to measure the spin of the Higgs boson. General helicity analyses prove the sensitivity of 
both decay modes to zero-spin of the Higgs boson, demonstrated by confronting J = zero-spin for illustration to 
J = 2 ± , i.e. spin = 2 and even/odd parity. 
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Research Foundation (NRF) funded by the Ministry of Education, Science and Technology (2012-0002746). The 
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